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Abstract 

It is demonstrated that, in the adiabatic approximation, non-Equatorial cir- 
cular orbits of particles in the Kerr metric (i.e. orbits of constant Boyer- 
Lindquist radius) remain circular under the influence of gravitational radi- 
ation reaction. A brief discussion is given of conditions for breakdown of 
adiabaticity and of whether slightly non-circular orbits are stable against the 
growth of eccentricity. 
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I. INTRODUCTION 



It has been shown some time ago that a particle in a circular orbit around a non-rotating 
black hole remains in a circular orbit under the influence of the gravitational radiation 
reaction arising from its orbital motion Although it has been suggested that the 
same holds true for "circular" orbits (meaning orbits of constant Boyer-Lindquist radius [Q) 
around rotating black holes, up to now it has not been shown to beyond first post-Newtonian 
order (see Ref. |^ for the post-Newtonian result, and see also Refs. and for more recent 
work) because of the difficulty of dealing with the little understood "Carter" constant of 
the motion, Q. To date no practical method has been developed for describing the rate 
of change of this constant due to gravitational radiation reaction for a generic orbit (see, 
however, Ref. 0), and without knowing this the evolution of the orbit cannot be predicted. 

In this paper, we study the relation between the rates of change of the three constants of 
the orbital motion (the energy, E, axial component of angular momentum, L, and the Carter 
constant, Q), for circular and almost-circular orbits. We first show that, for orbits which 
are precisely circular (in a sense which is well-defined within the adiabatic approximation), 
the rate of change of Q has just the value required for the orbit to evolve into a new circular 
orbit. Then, we extend the analysis to almost-circular orbits (to first order in the orbital 
eccentricity). This analysis leads to the result that in order for an initially-circular orbit 
to develop non-zero eccentricity, the back-reaction force (evaluated for the precisely-circular 
orbit) must resonate with the radial oscillations. Since in the case of a precisely-circular 
orbit the periodicity of the back-reaction force is determined solely by the ^-motion, we 
are led to the following conclusion. The only case in which an initially-circular orbit will 
develop an eccentricity is when there is a certain resonance between the (small-oscillation) 
radial motion and the 6'-motion. More specifically, this resonance condition is Tg = 2nTr , 
where n is an integer, and Tr and Tg are the (averaged) periods of the radial motion and 
the angular motion, correspondingly. 

Ryan |^ has recently examined circular orbits in the Kerr metric numerically and found 
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that the above resonance condition is never satisfied (for all black-hole and orbital parame- 
ters). We are thus led to the conclusion that orbits which are initially precisely circular will 
remain circular upon radiation-reaction evolution. We point out, however, that this conclu- 
sion does not address the issue of stability against the growth of a small initial eccentricity 
(this issue is further discussed in Sec. V below). 

This paper is organized as follows. In Section II we define instantaneous circularity of an 
orbit in terms of the instantaneous location and 4-velocity of the orbiting particle. We then 
derive a relation between the rates of change of E, L and Q for an instantaneously circular 
orbit that is perturbed by an arbitrary force and we show that this relation is precisely the 
one required for the circular orbit to evolve into a new circular orbit. One might naively 
interpret this result by itself as implying that an initially circular orbit will necessarily 
remain circular. However, in Section III we show that in order to predict the full evolution 
of initially circular orbits, it is necessary to carry the analysis (of the relation between 
the rates of change of the three constants of motion) to first order in the instantaneous 
eccentricity. This analysis, to first order in the eccentricity, is presented in Section IV, with 
the conclusion that initially circular orbits do, indeed, remain circular. Finally, in Section 
V we give some concluding remarks. 

II. INSTANTANEOUSLY CIRCULAR ORBITS 

In this section we shall define instantaneous circularity of an orbit and shall show that 
for any such orbits dQ/dr has just the right value so as to leave the evolving orbit circular. 
We shall show that this is true for any arbitrary force which acts on the orbiting particle (in 
fact, this result is precise, and is not limited to the adiabatic approximation). 

We take here the point of view that the radiative evolution may be viewed as the conse- 
quence of some "back-reaction force", which can be treated as any other external force [Q. 
We shall therefore begin by constructing a general formal expression for the evolution rate of 
all constants of (geodesic) motion, due to an arbitrary external force. Let C denote the con- 
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stant of motion in question. In Kerr, C may stand for either the energy E, the azimuthal 
angular momentum L, or the Carter constant Q (or any combination of these constants, 
hke e.g. the constant D defined below). We first express all these constants explicitly as 
functions of coordinates and covariant components of four- velocity 0, that is, 

C = C{x'^,u^). (1) 

For E and L, we simply take 

E = -ut , L = . (2) 

(Throughout, we use the standard Boyer-Lindquist coordinates, r,t,9,ip). The correspond- 
ing explicit expression for Q is slightly more complicated. We could use the familiar expres- 
sion based on the 6'-motion in Kerr: 

Q = u^ + cos^ 6 [a\l - ul) + sin^^ u^] . (3) 

We find it more convenient, however, to construct the expression for Q from the r-motion. 
It is straightforward to show that 

Q = [E{r^ + a^) - alf - {L - aEf - - Au^ (4) 

(this follows, for instance, from Equations (33.32b) and (33.33c) in Ref. fl^). Here, A = 
— 2Mr + a^, where M is the black hole's mass and aM is its angular momentum. For 
later convenience, we also write this equation in the form 

Q = H{r, E, L) - Aul , (5) 

where 

if (r, E, L) = A-^ [^(r^ + a^) - alf - {L - aEf - . (6) 

In view of Eq. (^, one readily sees that Eq. (Q) is just of the desired form, Q = Q{x'^,Ua) 
(with the simplification that the right-hand side does not depend on 6 or ug). 



When an external force is applied to the particle, C will evolve with time. To calculate 
its rate of change, we differentiate Eq. (|1]) with respect to the proper time r: 

^ = n^C, + E|^^- (7) 
dr a <^'^a dr 

Now, 

dUa 1 11 u /n\ 

^ = + l9,u,au''u'' (8) 

and 

DUr 



Dt 



Fa , (9) 



where D / Dt denotes covariant proper-time differentiation and is the force per unit rest 
mass. Equation (H) thus reads 



dC 
dr 



E^F.- (10) 



Notice that the term in brackets does not depend on the external force. When no external 
force is applied, C is conserved. Thus, the term in brackets must vanish identically. Equation 
(PUP therefore reads: 

f-^^^- (11) 

dr ^ dua 

This is the desired general expression for the evolution rate of all constants of motion. Note 
that this is the precise expression for the instantaneous rate of change of C. We have not 
used the adiabatic approximation (or any other approximation) so far. 

We now define an orbit to be instantaneously circular if its instantaneous values of E, 
L and Q [defined by Eqs. (0) and (||)] are precisely equal to those of some circular geodesic 
orbit. At a moment when the orbit (on which the force F^ acts) is instantaneously circular, 
Eq. (H) plus circularity {ur = 0) implies 

= -2Aur = . (12) 

OUr 



Inserting this into Eq. (0) gives 



dQ dQ dQ 

—r- = — -r* + -7: — r. 



dr dut * du, 



(13) 



Also, in view of Eq. (0), Eq. (|TT]) yields 



dE 



-Ft 



dL 
d^ 



(14) 



Substituting Eqs. (|) and ([l^) in Eq. (ITsI), we obtain 



^ dE ^ dL 

-f- = Q,E-\ h Q,L-r- ■ 

dr dr dr 



(15) 



Finally, using Eq. (^, we rewrite Eq. ([ISD as 




(16) 



This is our final expression for the actual, momentary, rate of change of Q due to the external 
force. [To avoid confusion , we emphasize that the partial derivatives in the right-hand side 
of Eq. ( P^ are to be calculated according to the explicit expression (^]. 

We come now to the second part of this calculation, that is to calculate the rate of change 
of Q (compared to that of E and L) required for taking a circular orbit into a new circular 
orbit. From Eq. (^ we obtain 



When applied to general geodesic orbits (with constant but arbitrary E, L and Q) this 
equation can be regarded as describing radial motion in an effective potential W{r, E, L, Q). 
Obviously that geodesic motion is circular if and only if the particle sits at a radius r where 
W = Q (or Ur = 0) and where Wr = (so the particle is at the minimum of the effective 
potential). Correspondingly, an orbit on which a force acts is instantaneously circular if 
and only if it instantaneously satisfies 



Aul = H{r, E,L)-Q = W^(r, E, L 



Q)- 



(17) 



W = Q 



(18) 
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We now let the force Fa continue to act, but only for an infinitesimal proper time, 5t. 
Following this action, the orbit will be characterized by new parameters, r', E' , L' and Q' . 
We denote the changes from the original parameters by 5, that is, 5r = r — r', 6E = E — E', 
etc. These changes are infinitesimal because Fa is allowed to act for only an infinitesimal 
time (5r, before the orbit is once more examined for circularity. The corresponding change 
in W is given by 

SW = W,r6r + W,eSE + W,lSL + W,q6Q . (19) 

In order for Eq. (18) to hold after the time 6t as well as before (i.e. in order for the orbit 
to remain circular), we must demand 

6W = 0. (20) 

We denote the value of Q which corresponds to a circular orbit (for given E and L) by 
Qcirc{E, L). Equations (p!9| ) and (^) thus yield, as a necessary condition for the orbit to 
remain circular 

W^rSr + W^eSE + Wx5L + W^qSQcuc = , (21) 
from which ^Qcirc is to be determined. Now Eq. ( P^D implies 

W^Q = ~1 , W,E = H^E , W^L = H^L, (22) 



which together with W^r = [Eq. ([TsD] reduces Eq. ( pTD to the form |11 



(JQcirc = H^eSE + H^l5L . (23) 

Finally, dividing by the infinitesimal proper time lapse 6t and taking the limit 6t 0, we 
obtain 

^ = He^ + H,^. (24) 

Comparing Eqs. (16) and (24), we arrive at the desired result. At any moment when 
the orbit, on which the arbitrary force Fa acts, is instantaneously circular, Fa produces an 
instantaneous evolution of the orbit's Cater constant given by 



dQ _ dQci 
dr dr 



circ 



(25) 



which maintains instantaneous circularity. 



For later convenience, let us define 



D = Q,^,-Q = D{E,L,Q) . 



(26) 



Circular orbits are thus characterized by = . Equation (^) then implies that if at a 



particular moment D = , then 



dDldr = . 



(27) 



It should be emphasized again that all the calculations done so far are precise, and do 
not depend on the adiabatic approximation. Note also that these calculations refer to the 
instantaneous rate of change of the constants of motion, at a moment when the orbit is 
instantaneously circular. As we shall see in the next section, Eq. ( ^Tf ) by itself does not 
imply that an initially-circular orbit will necessarily remain circular. One must go to the 
next order in the eccentricity in order to derive this result. 



Equation ( ^71) (which holds whenever D = 0, i.e. whenever the orbit is instantaneously 
circular) has a trivial exact solution, D{t) = . Does this necessarily mean that an in- 
stantaneously circular orbit will remain circular forever? We shall immediately see that, in 
principle, the answer is no (though, we shall also see later that, within the adiabatic ap- 
proximation, in most cases a circular orbit will remain circular). In fact, the trivial solution 
D{t) = to Eq. (|27| ) is physically meaningless. 

To illustrate this, consider the analogous (but much simpler) problem of a free particle 
in one-dimensional Newtonian mechanics. The particle's energy (per unit rest mass) is 
K = (1/2) [dx/drY . Assume now that a constant external force F (per unit rest mass) is 
applied on the particle. The evolution of K is then given by [in analogy with Eq. (pH]) ] 
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or, in terms of K itself, by 

dK 



V2FK^/^ . (29) 
dr 



Now, assume that the particle is initially at rest, i.e. K = . Then, from Eq. (p9D, 
dK/dr = . Equation (^) then admits a trivial exact solution, 

K{t) = . (30) 

This trivial solution, which means that the particle will remain at rest forever, is obviously 
wrong. Instead, the particle will certainly accelerate, and the true physical solution is [l^ 



K = (FV2)r2 . (31) 

The situation here with respect to Eq. (^) and its unphysical solution -D(r) = is just 
analogous. As will be shown below (Section IV), in the extension of Eq. (^T]) to slightly- 
eccentric orbits, dD/dr is (to the leading order) proportional to D^^"^ . Then, in addition 
to the trivial solution D{t) = 0, there exists a non-trivial solution in which for some short 
period D grows like [analogous to Eq. (|31|)], and this is the physical solution. This implies 
that, momentarily, the eccentricity (which is proportional to D^^'^) will grow linearly with 
r, even if initially it vanishes precisely. 

Recall, however, that we are not particularly interested here in the momentary rate of 
change at a specific point along the orbit. Rather, we are interested in the effective long-term 
evolution of D . Within the adiabatic approximation, this long-term evolution is described 
by an equation of the form 

D = S{E,L,D). (32) 

Hereafter, an overdot denotes the long-term rate of change (with respect to the proper- 
time, r), obtained from the momentary equation of motion by averaging over many periods 



(in Section IV we shall describe this averaging procedure in more detail). The long-term 
evolution of circular orbits will depend on the asymptotic behavior of S near D = 0. It is 
not difficult to show, based on Eq. (p7D, that the zero-order term [i.e. S{D = 0)] vanishes 
identically |T^. As we shall see in Section IV below, the general asymptotic form of S is 



given by 

S{E, L, D) = Si{E, L)D^/^ + S2{E, L)D + O^D"^/^) . (33) 

The evolution of an instantaneously circular orbit will depend crucially on whether Si van- 
ishes or not. In the case Si = , we can approximate Eq. (|33D by the linear equation 



D = S2D . (34) 

Then, an initial value D = ensures that D will remain zero forever (we are not concerned 
here about stability to small initial perturbations, though we shall discuss the question of 
stability briefly in Sec. V). On the other hand, if 5*1 is non-zero, we can approximate Eq. 
©by 

D = SiD^'^ . (35) 

As was explained above, the trivial solution D{t) = is physically meaningless. In that 
case, the physical solution will be (at least as long as 5*1 > 0) 

D{t) = (Si/2)V2 , (36) 

which means that the an instantaneously circular orbit will evolve into an eccentric one. 

The above considerations make it clear that the value of Si is crucial for our discussion. 
An orbit which is initially instantaneously circular will (in the adiabatic approximation) 
remain circular if and only if 5*1 = . In the next section we shall calculate 5*1, and show 
that it generically vanishes. Only orbits which satisfy a certain resonance condition may 
have non-zero Si . 
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IV. SLIGHTLY ECCENTRIC ORBITS 



We now analyse the adiabatic evolution of D for slightly eccentric orbits, to leading 
order in the eccentricity. In order to use Eq. (|TT|) for the calculation of dD/dr, we must first 
express D in the form D{x^,Ua)- To simplify the notation, we make use of Eq. (j^) , and 
simply write L instead of u^p and —E instead of Ut . Recalling that Qdrc = QcndE, L), we 
have from Eq. (^) 

D = QcUE, L) - Q{r, Ur, E, L) = D{r, Ur, E, L) . (37) 

Here, Q{r,Ur, E, L) is to be understood in terms of Eq. (1). This function is analytic in 
its arguments, and it is also not difficult to show that (except perhaps at the "last stable 
circular orbit", which must be excluded here), Qchc{E, L) is analytic in E and L. Therefore, 
D{r,Ur, E, L) in Eq. ( pTj) is also analytic in its arguments. Now, Eq. ( pA]) yields 

^ = -D^EFt + D^lF^ + D^^Fr = D'Ft + D^F^ + Fr = F, , (38) 

where the index j runs over the three coordinates r,t,(p , and 

= dD/duj . (39) 

Equation (|38D [like Eq. (p!T])j describes the precise, instantaneous rate of change of D, for 
any orbit on which any force Fa is acting. We have not made use of the adiabatic (or any 
other) approximation so far. 

From the analyticity of -D(r, Ur-, E, L) it follows that the functions 

= D^{r,Ur,E,L) (40) 

are analytic in their arguments. Now, from the validity of Eq. (^) for an external force 
of any type, it follows that for an instantaneously circular orbit [i.e. when Ur = and 
= '"circ(-£', L)], all three functions vanish. In other words, for any given E and L, 

D^{r = ro,Ur = 0,E,L) = , (41) 
11 



where tq = rarc{E,L) . We now expand D^{r,Ur, E, L) around (r = rQ,Ur = 0). In view of 
the analyticity of these functions, we have 

D^{r, Ur, E, L) = SrD^^^ + UrD%^ + 0{Sr\ ul, UrSr) , (42) 

where 6r = r — rQ , and the functions Z)-(, and D^^^ (which, again, are analytic functions of 
r, Ur, E, L) are evaluated at (5r = 0, Ur = 0) . 

We are now going to use two approximations (or expansions): the adiabatic approxima- 
tion, and the small-eccentricity approximation. These two approximations are unrelated, 
and should not be confused with each other. The exact instantaneous equation of motion 
of D is [cf. Eq. (|3|)] 

§^D'F,. (43) 

The adiabatic approximation means that the external force Fj is assumed to be small, and the 
right-hand side is to be evaluated to linear order in it. The small- eccentricity approximation 
means that the eccentricity is assumed small and Eq. ( ^31) is evaluated to first-order in it. 

In view of the small-eccentricity approximation, the higher-order terms in the right-hand 
side of Eq. (^) can be omitted. Substitution in Eq. ( ^31 ) then yields 

^ = [5riT)D^^r + UriT)D%^] F,{t) . (44) 

Let us now examine the implications of the two approximations used here on the expres- 
sion in the right-hand side. In view of the adiabatic approximation, the term in brackets 
is to be evaluated as if the constants of motion are fixed and the motion is geodesic. In 
the most general case, and D^^^ are (like D) functions of {r,Ur, E, L). Here, due to 
the adiabatic approximation, we can fix E and L . Moreover, since Sr and Ur are already 
first-order in the eccentricity, when evaluating D^^ and D^^^ we may take Ur = and r = tq 



. Thus, in Eq. (0), D^^. and are just constants (which depend parametrically on E and 
L) 0. 

Turn now to evaluate (5r(r) and Urir) in the right-hand side of Eq. (^). Like the entire 
term in brackets, they are to be evaluated as if the motion is geodesic (with fixed E, L, D). 
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In view of the small-eccentricity approximation, we need only calculate (5r(r) and 'Ur(r) to 
the leading order in the eccentricity. We start from the "effective-potential" relation 

Aul = W{r, E, L, Q) (45) 

[cf. Eq. (p!7|)], and, recalling that 5''"'' = g^^ = A/p^ , we write it as 

{dr/drf = {^VpVr = P~^W/\ . (46) 

Here, 

p2 ^r^ + a^cos^^. (47) 



In Eq. (^61) , as it stands, the radial motion is coupled to the ^-motion, through p . In order 
to decouple the two motions, we define a new independent variable A by 

dX/dT = p-^ . (48) 

The radial equation of motion now becomes 

{dr/dX)^ = WA , (49) 

in which the right-hand side is independent of 6 . 



From Eqs. (]T7|) and (|g), W is an analytic function of {r, E, L,Q) . Writing Q = 
Qcirc{E,L) — D , and recalling the analyticity of Qcirc{E,L) , W can be expressed as an 
analytic function of (r, E, L, D) . For an instantaneously circular orbit (i.e. for r = tq and 
D = 0), both W and W^r vanish [cf Eq. ([T8|)j. The expansion of W near instantaneous 
circularity must therefore be of the form 

W = Ad + Bdr"^ + 0{D\D6r,6r^) , (50) 

where A and ]3 are some functions of E and L . Correspondingly, the expansion of the 
right-hand side of Eq. (HDI) will take the form 



WA = AD - B6r^ + 0{D^,D6r,6r^) , (51) 
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where A = AAq and B = — -BAq , and where 

Ao = A{r = ro{E,L)) . (52) 



Combining Eqs. (|49|) and (^T]), and restricting attention to the leading-order eccentricity 
effect, we obtain 

{d5r/d\f = AD- B5r^ . (53) 

Note that A and B are some functions of E and L only. Therefore, as explained above, they 
can be regarded here as fixed parameters. 

Equation (^) describes a simple harmonic oscillator. Its general solution is 

(5r(A) = kVd cos[w,(A - Aq)] (54) 



where K = ^ A/ B and oj^- = \fB are parameters that depend on E and L only. (Do not 
confuse K here with K of Sec. III). Using this result to calculate m^, we find 

Ur = Qrr-T = ^ ""^"TT ^ — K y/DuJr Sm[Ur{X — Aq)] . (55) 

dr dX 

To simplify the notation, we shall hereafter absorb the constant Aq into A (by shifting the 
origin of the latter if necessary). Substituting Eqs. (Q) and ( p^D into Eq. (^Ij), we obtain 
dD 



dr 



KVD (D;'^ COs(cJ^A) - D%^A-^UrSm{Ur\))] Ej{T) . (56) 

From Eq. (|56|) it is already clear that, so far as the instantaneous evolution is concerned, 
dD/dr is indeed proportional to \/D . Thus, as explained in Section III, although D = Q 
yields dD/dr = 0, an instantaneously circular orbit will not remain circular later on. Instead, 
a momentary growth of D like is to be anticipated. However, we are primarily interested 
here in the long-term adiabatic evolution of D . The latter is to be obtained from Eq. ( |56D 
by averaging it in r over many periods. In order to perform this long-term averaging, we 
must first take a closer look at the nature of the time dependence of the force Fj . 

Since the term in brackets on the right-hand side of Eq. (^) is already proportional to 
(i.e. to the eccentricity), when evaluating Fj(T) we are allowed to assume that the orbit 
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is a precisely circular geodesic. The 9 motion of such an orbit is periodic in r. Therefore, 
the backreaction force must be periodic also. The various points along the circular geodesic 
orbit are physically distinguishable from one another only by the values of 6 and dO/dr. It 
therefore follows that after completing a full cycle of the ^-motion, Fjij) will return to its 
original value. A closer look, however, reveals that, because of the reflection symmetry of 
the Kerr geometry, the ^-motion is symmetric with respect to the equatorial plane. As a 
consequence, the period of Fj{T) will in fact be half of that of the full ^-motion cycle. 

In order to facilitate the calculations, it is convenient to transform Eq. ( ^6|) from r to A. 
Using Eq. (^8]), we obtain 

dD 



(z^f^cosKA) - D%^A~^UrSm{uJr\))] Fj{T{\)) . (57) 



According to our expansion scheme, we need only evaluate the term in brackets to zero order 
in the eccentricity. That is, we can replace p and A by their circular counterparts, po and 
Aq , where Aq is the constant defined in Eq. (|5), and po is a function of 9 only, defined by 

pl = rl + cos^ 9 . (58) 

We therefore obtain 

^ = K\fD [pi, cos(c^,A) - /^l, Ao 'a;, sin(a;,A)) [p^F,] . (59) 

The term in brackets depends on r (and hence on A) through its dependance on the 9- 
motion. It is obvious from Eq. ([58| ) that po is periodic in r, again with a period which 
is just half that of the 6'-cycle. Consequently, the entire term in brackets is also periodic 
(with that one-half 9 period). Let us examine now the periodicity of this term with respect 
to A. Again, the ^-motion is periodic in A , and the reflection symmetry implies that the 
motion at < 7r/2 is just symmetric to that at 9 > 7t/2. (This can also be deduced directly 
from the fact that, in the equation of motion for 9{X), {d9/dXf' = 0(6*) (cf. Ref. ||10[)) the 
function <d{9) admits a reflection symmetry about 7r/2). Thus, if we denote the A-period of 
the ^-motion by 
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Ke = 271 /cue , (60) 

then the term in brackets in Eq. (^) has a A-period of A0/2 . Correspondingly, the Fourier 
transform of this term will take the form 



pIFj = J2 Cy^"^'^ + C.C. , (61) 

n=0 



where C.C. means the complex conjugate of the preceding term. 
Substituting this expansion in Eq. ( |59| ) yields 

dD 



00 



D [K [D], cos(cj,A) - Df„^Ao ^cu, sin(cu,A))J ^ G^^e^"*"*^ + C.C. . (62) 

n=0 



It is convenient to transform the last expression from Sine and Cosine to exponential func- 
tions. The term in brackets can be expressed as 

where 

Ki ^ {K/2) [d], ± iD^^^^^-,'uJr) . (64) 
Eq. ( |62D then takes the form 

— = Vd (i^ie*^--^ + Kie"*"'-^) J2 G"e2"*"«^ + C.C. . (65) 

n=0 

Defining now 

Gl = KiC] , (66) 

we obtain 

dD °° 

— = ^/DY, (G'!^e^(2'^'^»+^')^ + G^e^^^""^*"^')^) + C.C. . (67) 
^-^ n=0 

Recall that in this equation the coefficients G^ depend on E and L , but not on A . 

Equation ( pTD describes (within the adiabatic approximation, and to leading order in the 
eccentricity) the instantaneous rate of change of D . In order to obtain from it the long-term 
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rate of change, we simply need to take its average over a sufficiently long period of A . To 
that end, for any function U (A) , we formally define the long-term averaged rate of change 

/dU\ AU 

( — ) = lim — , (68) 

\dX / AA^oo AA ^ ^ 

where AU and AA denote the difference in U and A , correspondingly, between the two 
extremes of the A-interval considered. Although the averaging is over times long compared 
to l/uj0 and 1/uJr, it is still short compared to the radiation reaction timescale, which 
is the time for substantial orbital inspiral. (Recall that since we are using the adiabatic 
approximation here, if U also depends on the "constants of motion" , they must be taken as 
fixed constants in this averaging process.) The averaging of the right-hand side of Eq. (|67D 
is trivial, in that the term e«{2n<^e±'^r)A ^jjj average to zero, unless 2nujg ± cj^ = , in which 
case it averages to one. Since both ug and Ur are taken to be positive, we need only worry 
about those terms with 2nuJe — Ur in the exponent. We thus obtain 

(^ — j = ^Y.Gn5i2n.,-..), (69) 

where 5 is a function whose value is unity when 2nujg — Ur = and zero otherwise, and 
Gn = 2^ (^-)' ^ meaning the "real part of". 

At this stage it is already clear that, unless a certain resonance condition is satisfied 
{uJr = 2nuj0 for some ra), the right-hand side of Eq. ( p9| ) will vanish. Before we further 
discuss the meaning and implications of this resonance condition, however, we shall more 
directly connect our result ( |69D to the notation used in Section III, and in particular to the 
parameter Si . Equation ( |55| ) is to be obtained from the momentary rate of change of D 
by averaging over proper time. For any function U{t)^ the long-term proper-time average 
(denoted by an over dot) may be formally defined as 

U = hm -— = hm — lim — — = J M — ) . (70) 

At^oo At aa^oo At aa^oo AA \ dX 



where 
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is a constant that depends on the orbit. Applying this procedure to we obtain |T5 



/ fin\ °° 
D = l^ — j = ^/DJ-^ YlGJi2n.,-u,.) . (72) 

Comparing now Eq. ( [7^ ) to Eq. (^), we find 

oo 

Si = J Gn5(2nujg-UJr) ■ (73) 

n=0 

The imphcations of this result for the long-term evolution of D are obvious. There are two 
different cases: 

a) The resonant case: there exists an integer n such that uo^ = 2nujg . In that case, we 
have 

5-1 = J-'Gn , (74) 

which is likely to be nonzero in the general case. Then, as discussed in Section III [cf Eq. 
(p6|)] , D will grow like r^, which means that the eccentricity will grow linearly with r. 

b) The non-resonant case: there exists no integer n for which Ur = 2nujg . In that case, 
we simply have 

^1 = 0, (75) 

and the equation of evolution will read 

D = S2D + 0{D^/^) . (76) 

[cf Eq. (|33|) ]. In this case an orbit which is initially precisely circular will remain circular 
(within the adiabatic limit). (The sign of 5*2 will determine the stability against growth of 
small initial eccentricity). 

For resonant orbits, we have 

Ae = 2nK (77) 

where Kg and are the A-periods of the 9- and r-motions, correspondingly. It would 
sometimes be convenient to translate this expression to t-periods. One finds that, not 
surprisingly, the resonance condition is 



Te = 2nTr , (78) 

where Tg is the t-period of the ^-motion, and is the averaged t-period of the radial motions. 
[The radial motion, expressed in terms of t (or r), is quasi-periodic rather than periodic, 
because it is modulated by the 6'-motion. The 6'-motion itself is periodic in either t or r — 
first, because we are considering a circular orbit here, and second, because the resonance 
condition (^) ensures that each time 9 returns to its original value, r does also (but not 
vice versa).] 

V. CONCLUSION 

We have shown that, within the adiabatic approximation, an orbit which is initially 
precisely circular will remain circular, under the action of the radiation-reaction force. The 
only exception is if the orbit satisfies the resonance condition Tg = 2nTr , for some integer 
n, where Tg is the 6'-motion period and is the (averaged) period of the small-oscillation 
radial motion. However, circular orbits around Kerr never satisfy this resonance condition 
We therefore conclude that, within the adiabatic approximation, an orbit which is 
initially circular will remain circular. 

There are two caveats which should be mentioned here. First, no attempt has been 
made to address the issue of stability against the growth of a small initial eccentricity. This 
stability would depend on the sign of the coefficient 5*2 in Eq. (|3^) above, which was not 
calculated here. 

Second, our conclusion that circular orbits must remain circular is only valid within the 
adiabatic limit, i.e. in the limit /i/M — 0, where /i is the mass of the small object. In reality, 
since the ratio fi/M is always finite, an initially circular orbit will develop some eccentricity. 
For concreteness, consider a particle with fi <^ M which at some initial stage (which we 
denote stage 1) moves along a circular orbit with Boyer-Lindquist radius ri. Later on, the 
orbit shrinks due to radiation reaction, until (at stage 2) it passes through a new radius, 
r2 < Vi. Then, for every finite fJ^/M, one should expect non-zero eccentricity to be present at 
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stage 2. The above analysis, however, imphes that the eccentricity at stage 2 will decrease 
with yu (for fixed ri, r2 and M), and will vanish at the limit n/M —>■ (presumably like 

The small eccentricity of non-adiabatic origin mentioned above could in principle seed 
an exponential growth of eccentricity if 5*2 > 0. In such a situation, an initially circular 
orbit may evolve into a very non-circular one (even for small fi/M). The feasibility of this 
scenario depends, of course, on the relevant values of fi/M and on whether 5*2 > and also 
on the available range of r- values (over which 5*2 is positive). 

In reality, we know that in the Newtonian limit the orbits become more and more circular 
as they shrink due to radiation reaction, so we expect ^2 to be negative as long as r ^ M 
(recall that the Newtonian approximation should hold at r ^ M even if the black hole is 
spinning). Consequently, the range over which 5*2 might be positive is bounded. We can 
therefore expect that if /i/M is sufficiently small (and if r2 is not too close to the "last stable 
circular orbit" fl^), the instability will not have enough time to build up, and an initially- 
circular orbit will indeed remain circular throughout the inspiral, to a good approximation. 
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